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Abstract 

We consider a classical lattice gas that consists of more than one "species" of parti- 
cles (like a spin-| Ising model or the atomic limit of the extended Hubbard model), 
whose ground-state phase diagram is macroscopically degenerate. This gas is cou- 
pled component- wise and in the Falicov-Kimball-like manner to a multicomponent 
free-fermion gas. We show rigorously that a component-wise coupling of the classical 
subsystem to the quantum one orders the classical subsystem so that the macro- 
scopic degeneracy is removed. 
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1 Introduction 

The Hamiltonians of the so-called strongly correlated electron systems de- 
scribe in fact multicomponent gases of quasiparticle excitations in a crystal. 
The components are labeled, by orbital and/or spin degrees of freedom, for 
instance. As a rule, revealing cooperative phenomena in such systems, like 
long-range orders, is a hard task. Then, in the course of studies, one simpli- 
fies the Hamiltonian by discarding some labels of quasiparticles, setting some 
coupling constants to zero and sending other to infinity. Simultaneously, one 
looses possibilities of accounting for some physics. But once a sufficient level 
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of understanding of a simplified model is reached, a reverse process starts. 
This might be a brief history of the Falicov-Kimball model (and only to some 
extent of the Hubbard model) studies. 

The original model put forward by Falicov and Kimball [1,2] refers to two 
groups of quasiparticles excitations, called localized holes and itinerant elec- 
trons, with the members of each group labeled by band and spin indices. 
Only a few, out of many, on-site interactions of these quasiparticles are re- 
tained and the strength of the on-site hole-hole repulsion is sent to infinity 
[2]. The extreme simplification of this model, that consists of one band of spin- 
less electrons and of spinless holes, is known as the spinless Falicov-Kimball 
model. In contradistinction to the case of Hubbard model, it has been possi- 
ble to demonstrate, without uncontrollable approximations, that the spinless 
Falicov-Kimball model exhibits a staggered long-range order [3,4,5]. Subse- 
quent research extended considerably this result and supplied us with a fairly 
broad knowledge of the phase diagram. Naturally, the strongest results have 
been obtained at zero temperature (see [6,7] for more details and an extensive 
list of references). 

Then, one observes a resurgence of interest in effects of additional degrees of 
freedom and additional interactions between localized particles. Most atten- 
tion has been paid to effects related to spin degrees of freedom. Thus models, 
where the localized and intinerant quasiparticles carry spin-| [8,9,10,11,12,13] 
have been investigated. In [8,9,13] a finite repulsion between localized particles 
with opposite spins has been included. Quite recently, it was argued that an in- 
teresting physics could be accounted for by multicomponent Falicov-Kimball 
models [14]. As a preliminary step, a kind of Falicov-Kimball model with 
2s + 1 components of itinerant particles and 25* + 1 components of localized 
ones was solved exactly in the limit of infinite dimensions. Concerning inter- 
actions between localized particles, only an on-site repulsion between different 
components of localized particles was included and its strength was sent to 
infinity, thereby restricting the space of states available for the system. Unlike 
in the quoted above papers on the spin-| case, s may not be equal to S. 

The multicomponent Falicov-Kimball models, we referred to above, belong to 
a class of systems whose generic Hamiltonian is of the form 

Ha = E T ° + E ^X,X,, + H i- (!) 

The first sum, over cr-components of itinerant fermions (later on called briefly 
electrons), accounts for their hopping. The second sum, over lattice sites x, 
cr-components of electrons and ^-components of localized particles (later on 
called ions), represents the usual Falicov-Kimball-like coupling between the 
itinerant and localized particles. With n e x a , n l xr) being the occupation number 
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operators at site x of a-electrons and 77-ions, respectively, it contributes the 
energy U ar) if and only if the site x is occupied by a cr-electron and an rj- 
ion. The third term, Hi, stands for an interaction between ions, diagonal in 
their occupation-number representation. Therefore, the system described by 
the Hamiltonian H\ can be looked upon as a multicomponent classical gas (of 
ions) coupled to a multicomponent free-fermion gas (of electrons). 

In the main body of this paper we study an instance of (1) with a component- 
wise coupling between the electrons and ions, i.e., we set U arj = U a 5 arj . The 
component-wise coupling may appear naturally in situations, where the com- 
ponents differ significantly in their physical characteristics. Actually, we make 
the coupling component-independent, U a = U for all a , and assume equal 
number of electrons' and ions' components, for simplicity. If Hi involves only 
finite-range interactions that admit a representation in terms of a m-potential 
[15], then some aspects of the ground state of the model, defined by H A , con- 
cerned with the distribution of ions on the underlying lattice, can readily be 
investigated by the methods developed in [4,16]. 

However, if the itinerant and localized particles are interpreted as different 
sorts of electrons, as is the case in the original paper by Falicov and Kimball, 
and the coupling between electrons stands for a screened Coulomb repulsion, 
then an ion-independent coupling, i.e., U ar] = U a , becomes more appropriate. 
Such a coupling has been studied rigorously in [8,9] and in the limit of infinite 
dimensions in [14]. The multicomponent Falicov-Kimball models with ion- 
independent coupling will be the subject of our subsequent paper. 

In the next section we provide necessary details and properties of our model 
with the component-wise coupling, as well as definitions useful in considera- 
tions of ground-state phase diagrams in the grand-canonical ensemble. Then, 
in Section 3 we study the ground-state phase diagram at the hole-particle 
symmetry point in the plane of chemical potentials of electrons and ions, for 
arbitrary strength U of the electron-ion coupling. After that, in Section 4 we 
go beyond the hole-particle symmetry point, but only in the strong-coupling 
regime. Finally, in Section 5, we provide a summary and discussion of the 
obtained results. 



2 The model and the grand canonical ground states 

In the sequel we study a version of the model (1) specified as follows. The 
underlying finite lattice A, consisting of | A | sites, is assumed to be a piece 
of a rf-dimensional hypercubic lattice Z d , such that it can be partitioned into 
two congruent sublattices: the even sublattice, A e , and the odd sublattice, A . 
Each site of A has z = 2d nearest neighbours. Then, the nearest-neighbour 
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sites of a site in A e belong to A° and vice versa. We shall use the notation 
(x, y) = (x, y)i, (x, y)i, % — 1, 2, . . . , x, y e A, to denote the first and the i-th 
nearest neighbours on A, that is the unordered pairs of sites whose Euclidean 
distance in the units of the lattice constant are l,-s/2,2,... , for i = 1,2,3,..., 
respectively. Consequently, Y^i x ,y)i w * n stand for the summation over all the 
i-th nearest-neighbour pairs of sites, with each pair counted once. 

The quantum subsystem, consisting of r ^ 2 different kinds of itinerant par- 
ticles (electrons) labelled by a — 1, 2, . . . , r, is described in terms of creation 
and annihilation operators of a a-electron at site x G A: c+ CT , c Xj(7 , respec- 

tively. The operators \^c^ a ,c x ^j ' ' satisfy the canonical anticommutation 
relations. Then, the hopping of a-electrons is specified by: 

^ = *E (4/^ + h-c) , (2) 



with t being the nearest-neighbour hopping intensity. Another relevant observ- 
able is the total electron-number operator N e = X^A^, where a-electron- 
number operators iV e)<7 = J2 X n t,a an d n x ,<r = c x,a c x,a ■ 

The classical subsystem consists also of r ^ 2 different kinds of particles (lo- 
calized ions). The states of this subsystem can be described by a collection of 

r . *| a=l,...,r 

particle-occupation numbers C = \ n l x ^ j eA , with n l X(T = 0, 1 (no two ions 
occupying the same site can be of the same kind), called the ion configura- 
tions. The total number of ions is iVj = J2a ^i,<r> where the numbers of a-ions 
Ni t a = J2x n xa- I n contradistinction to the electron subsystem, we introduce 
a direct ion-ion interaction. Each pair of ions occupying a site contributes the 
energy 21. Moreover, two ions that occupy two nearest-neighbour sites repell 
each other, contributing the energy W/2 > 0. Thus, the ion-ion interaction 
Hamiltonian, Hi, reads 

^ = wEE + t £ « ( 3 ) 



where n l x = J2a n x ,a an d Ha'<a" stands for the summation over unordered 
pairs of indices (a', a"), each pair being counted once. 

Clearly, in the composite system, whose Hamiltonian is given by (1), (2) and 
(3) (with arbitrary matrix U arj ), the particle-number operators iV e)(T (hence N e ) 
are conserved. The peculiarity of this system is that the site-occupation num- 
bers n x a (hence N^ a and Nj) are also conserved and therefore the description 

( ■ 1 <r=l,...,r 

of the classical subsystem in terms of the ion configurations C = \n xa J ^ 
remains valid. 
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Finally, the coupling between the classical and quantum subsystems is chosen 
to be the component-wise one, with 

U av = 2U5 av . (4) 



This concludes the specification of our system. 

In what follows, we shall study the ground-state phase diagram of the sys- 
tem defined by (l)-(4), in the grand canonical ensemble. Ultimately, this task 
amounts to determining at each point (/2 e , of the plane of chemical poten- 
tials of electrons and ions, respectively, the lowest energy eigensubspace of the 
Hamiltonian 

Ha (fie, P-i) = H A - fx e N e - /^A^. (5) 



Actually, our purposes will be more restricted. Namely, we shall be interested 
only in determining the set of periodic ground-state configurations of ions at 
each point in the (// e ,/2j) plane. 

In studies of grand-canonical phase diagrams an important role is played by 
unitary transformations that exchange particles and holes: rf xa — > 1 — rf xa 
and n l xa — > 1 — n X(7 , and for some leave the Hamiltonian Ha (jj> e , P>i) 

invariant. For the electron subsystem of (5) such a role is played by the trans- 
formations: c xa — > (— l)' x 'c XjCr , with \x\ = \(xi, . . . , Xd)\ = J2i=i\ x i\- At the 
hole-particle symmetry point, the system under consideration has 

very special properties, which facilitate studies of its phase diagram (among 
them is the property that the grand-canonical expectations of n X(7 ,n X(7 as- 
sume the value of |, independently of temperature). One easily finds that this 
system is hole-particle invariant at the point (ji®,/!®), where 

£ = U, tf = is° + U, %=™rz + I(r-l). 



For future purposes, it is convenient to introduce the shifted chemical poten- 
tials 

fi e = fi e - jl e , Hi = iii- fit, 



the new Hamiltonian i?A (/-te, /-0 = Ha (p> e , p>i), and express i?A (/x e , /ij) as 
follows 

H A (/i e , = £ (H™ - fi e N^ - frNi,,,) + Hi (6) 



5 



x 

= T ff + 2tf£(n^-±) (7) 
H^Hi-vfNi 




The second equalities in (7, 8) hold up to a term which is independent of 
the particle states. The Hamiltonians H% K and Hf are invariant with respect 
to the hole-particle transformation, consequently H\(fa,fa) is hole-particle 
invariant at the point fa — fa — 0. 

We define Eq (fa, fa) to be the ground-state energy of Ha (fa,, fa) for a given 
configuration of ions C. Then the ground-state energy of the Hamiltonian 
H A (fa,fa), denoted E G (fa,fa), is Eg (fa, fa) = mm{E c (fa,fa) : C}, and 
the minimum is attained at the set G of the ground-state configurations of 
ions. 



3 Ground-state phase diagram at the hole-particle symmetry point 

In this section we shall describe the ground-state configurations of ions in the 
systems given by Hf and H\ (0, 0). 

Let us consider first the classical subsystem of ions, given by Hf. The nature 
of its ground-state configurations of ions is explained by the Theorem 1, which 
follows. Briefly, there are two domains of couplings /, W, differentiated by the 
ratio 5 = AI/zW. For 5 < 1, there are exactly two ground-state configurations, 
while for 5^1 there is a macroscopic degeneracy (i.e., the number of ground- 
state configurations grows exponentially with the size |A| of the underlying 
lattice). 

Theorem 1 Consider the system given by Hf , with W > 0. 

(i) If 4/ < zW , the set of ground-state configurations consists of exactly two 
ion configurations with the following property: each site of one of the sublat- 
tices of A is occupied by r ions while the complementary sublattice is empty. 

(ii) If 47 > zW , the set of ground- state configurations contains macroscopi- 
cally many elements. 

For even r,r — 2k, any ion configuration with k ions in each site is the 
ground-state one; the degeneracy of the ground state is (Cy)' A ', where = 
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r\/k\{r — k)\. 

For odd r,r = 2k+l, any ion configuration satisfying the condition: each site 
of one sublattice is occupied by k ions while each site of the complementary 
sublattice by (k + 1) ions, is the ground-state configuration; the degeneracy is 

(c r fc )i A i/ 2 (cf +1 ))i A i/ 2 . 

(Hi) If 41 = zW the set of ground-state configurations contains macroscopi- 
cally many elements: any configuration of kind (i) or (ii) is the ground-state 
configuration. 

A proof of Theorem 1 is given in appendix. 

Now, consider the composite system that consists of r components of clas- 
sical ions coupled to r components of quantum electrons, described by the 
Hamiltonian H A (0,0). The following theorem characterizes the ground-state 
configurations of ions at the hole-particle symmetry point. 

Theorem 2 Consider the system given by Ha(0,0) with U ^ and W > 0. 

(i) If 41 < zW , then the ground-state energy is attained for exactly two ion 
configurations that satisfy the following condition. The ions fill out completely 
one sublattice of A, i.e., each site of this sublattice is occupied by r ions, while 
the complementary sublattice is empty. 

(ii) If 41 > zW , then the set of ground-state configurations of ions contains 
finitely many elements; any two nearest-neighbour sites cannot be occupied by 
ions of the same kind. 

Moreover, for r even, r = 2k, there are exactly ground-state configurations 
of ions; each site of A is occupied by k ions. 

For r odd, r = 2k + 1, there are exactly 2C^ ground-state configurations of 
ions with the following property: each site of one sublattice is occupied by k 
ions, while each site of the complementary sublattice contains k + 1 ions. 
(Hi) If 41 = zW , the set of ground state configurations of ions consists of 
finitely many elements: any configuration of kind (i) or (ii) is the ground- 
state configuration. 

Thus, the effect of the component- wise coupling of our classical gas of ions with 
quantum electrons amounts to removing the macroscopic degeneracy of the 
domain of /, W, where 5^1. The domain, where 5 < 1, of finite degeneracy 
is left unchanged. This effect does not depend neither on the hopping intensity 
t nor on the strength U of the ion-electron coupling. 

To prove the Theorem 2, let us note that since the operators H^ K (a = 
1, . . . , r), H®, commute pairwise, for given configuration C 

E c (0,0) = j:E™ + Ei J , (9) 

a 

where is the ground-state energy of H^ K for the configuration C a = 
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{^x,cr} gA that is the restriction of C to cr-ions, and E l c is the value of if° at 
C. Let Eq K be the ground-state energy of H™, i.e. = min {Eg? : C CT } 

and the minimum is attatined at the configurations C a G G„ K . The energy 
Eq K and the set G% K do not depend on a. Let = min {E l c : C} and let 
G l be the set of configurations for which the minimum is attained. In general 

£ (o,o)>5>£* + i&. (io) 

However, we shall show that the lower bound (10) of Eq K is attained for the 
set G of the configurations described in Theorem 2. 

Consider first the term J2a of (9). Clearly, 



min 




G 



and is attained for configurations in 

G FK = {(C 1 ,...,C r ) :C a eG™,a = l,...,r}. 



In order to determine G FK (and Eq K ) we follow [4]. Note that, for a fixed 
configuration C the operator H FK + UJ2 X ( n l,a ~~ i) i s the second quantized 
form of the one-particle operator h a = T + US a , where T and are |A| x |A| 
matrices with the following matrix elements [T] xy , [S a ] xy : 

{t, if x, y are nearest neighbours on A 
0, otherwise 



The lowest eigenvalue of H FK + U J2 X \ n x,a ~ i) i s the Fermi-sea energy of 
h(j corresponding to the Fermi level \x e = 0. Thus, in terms of the eigenvalues 
AJ, j — 1, . . . , |A|, of the operator h a 



E 



FK 



AJ<0 Z x z 



2 ^,<t 



-i [Tr|^| -Tr^] - | £ s x , ff = -^Tr|/i CT |, 
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since TrT = 0. Kennedy and Lieb have proved in [4] that 
-Tr\h a \ ^ -Tr (T 2 + U 2 ) 

and the lower bound is attained for exactly two configurations C a with the 
following property: one sublattice of A is completely occupied by <r-ions while 
the complementary sublattice is empty. Thus, Eq K = — |Tr (T 2 + £/ 2 ) and 
the described above configurations constitute G FK . 

The set G 1 is specified in Theorem 1. Now, since G FK n G % is nonvoid, the 
set G of ground-state configurations of H\(0 : 0) is G = G FK fl G 1 and the 
ground-state energy equals to E G (0, 0) = -§Tr (T 2 + U 2 ) + E l G . 

In the following section we shall show that the effect of removing macroscopic 
degeneracies of classical-ion phases persists off the hole-particle symmetry 
point. 



4 Ground states off the symmetry point in the strong-coupling 
regime 

Let us consider the system given by (6), (7) and (8) with ^ ^ 0. In order to 
simplify the arguments that follow, we set r = 2, with the label a interpreted 
as spin- 1 and taking values in {til}- Moreover, the calculations are carried 
out for d = 2, i.e. A is a piece of the square lattice, but the obtained results can 
be generalized to higher-dimensional hypercubic lattices in a straightforward 
way. 

As in the previous section, we start with the classical subsystem whose Hamil- 
tonian is — faNi, and is specified by (8). The ground-state phase diagram 
in the plane (I/W, pt-i/W) can be obtained by extending the calculations of 
Appendix A, and is shown in Fig. 1. 

The energies of ion configurations depend only on the occupation numbers 
of sites: n % x = 2,1,0. Consequently, the ground-state configurations in the 
domains B, C + , C_ are macroscopically degenerate due to the spin degeneracy. 

According to the Theorem 2, at the line /Xj = and in the domain B, the 
coupling to the electron subsystem orders the ion subsystem antiferromagnet- 
ically: each sublattice is occupied exclusively by one kind of ions. 

In order to determine the effect on the macroscopically degenerate ground- 
state configurations of the coupling to the electron subsystem, off the line 
\ii = 0, we resort to the strong-coupling limit. In this limit, the two components 
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Fig. 1. Ground-state phase diagram of H® — /xTVj, for fi ^ 0. The ground-state 
configurations inside the regions A, B, C+, D + are obtained by periodizing the 
elementary-plaquette configurations displayed inside the regions A, D + and those 
located by the starting points of the arrows in B, C+. The elementary-plaquette 
configurations located by the end points of the arrows in B, C + show the effect of 
the coupling to the electron subsystem. The second half of the diagram, for \i < 0, 
containing the second halves of the regions A, B and the regions C_, D_, can be 
obtained by reflecting the part with \i > in the line fi = and simultaneously re- 
placing the elementary-plaquette configurations, by their hole-particle counterparts 



■TV, 



n" 



— /i e A^ e o - of H\ (/i e , Hi) can be expanded into a power series in 1/\U\ [16]. 
The obtained series is convergent absolutely and uniformly in |A|, provided 
\U\ > zt and \fi e \ < \U\ — zt (here z — 4) [17,18]. Due to symmetries of 
if a (fie, Hi) we can restrict our considerations to the case U > 0, and set 
N e , a = |A| - [4,16]. 



Let (/i e , fj,i), k — 1, 2, . . ., stand for the partial sums of the \j\J expansion 
of if a (/i e ,/ij). Then, up to a term independent of the ion configurations, the 
partial sum (fie, Hi) reads [16]: 



) (He, Hi) = -yE ( S *,T + s *d) + \ E S x,1 S x,l 

" X " X 

w _ , 
' t 2 9t 4 \ 
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+ 



' {x,y)i 



3t 4 t 4 
+ YQIP S ( S ^T S ?/,T + s *d%i) + ^3 2 ( s ^T s j/,T + 
(z,y)2 (x,j/) 3 
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Since H ( A ' (/i e ,/ij) differs from Hf — /ijiVj by a constant, at the zero order 
the phase diagram coincides with that of the classical subsystem. In the first 
order, (/z e ,/ij) is augmented by the Ising antiferromagnetic interaction of 
the strength t 2 /4U. Consequently, in the phase diagram of (// e ,//j) the 
domains B, C+, C_ are replaced by some domains C+\ whose 

boundaries are shifted with respect to the boundaries of B, C+, C- by terms 
of the order l/U. The set of ground-state configurations inside B^> consists 
only of the two antiferromagnetic configurations, while the sets of ground- 
state configurations inside C+\ remain the same as those in C + , C_. 
Passing to the next order, one finds that the phase diagram of H A {/i e , //j) 
contains some domains B^ 3 \ C+\ whose boundaries are shifted with 
respect to the boundaries of domains B^ l \ C+\ ci 1 ^ by terms of the order 
1/U 3 . In this order, the set of ground-state configurations inside B^ remains 
the same as that in On the other hand, the macroscopically degenerate 
sets of ground-state configurations in (Ci 1 ^) are replaced by exactly four 
configurations in C + (C_ ): one sublattice is doubly occupied (empty) while 
the complementary sublattice is ordered antiferromagnetically. Applying the 
arguments due to Kennedy [17] and Gruber et al. [18], we conclude that there 
are nonvoid domains B (co \ , such that C 5 (3) , Cf° } C C| 3) , 

q(p°) c £r(3)^ - ms ^ e ^ich the sets of ground-state configurations of H A (/i e , /ij) 
coincide with those inside B^ , C+\ respectively. 



5 Discussion 

The two-component classical system considered in Section 4 and given by the 
Hamiltonian Hf — is well known in the literature of strongly correlated 
electron systems as the atomic limit of the extended Hubbard model. A lot 
of effort has been devoted to studies of its phase diagram, see for instance 
[19,20] and references quoted there. When the two kinds of ions involved (in 
the context of Hubbard-like models they are interpreted as localized electrons 
with spin up or down) are allowed to hop, we arrive at the extended Hubbard 
model. This model is, since many years, a subject of vigorous research. One of 
the intriguing questions asked is what is the effect of hopping on the phases 
of localized electrons (determined in the atomic limit). It is widely known 
that the second-order perturbation theory, with respect to hopping, predicts 
that the ground-state configurations with only singly occupied sites (inside 
domain B) are ordered antiferromagnetically by a weak hopping. Brandt and 
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Stolze have shown [21] that the prediction for the phase with density | (|), 
inside domain C + (C_), is that it is ordered ferro magnetically by a weak 
hopping, what results in coexistence of charge and magnetic orders. Thus, 
according to the second-order perturbation theory the quantum fluctuations, 
introduced by the hopping, remove the spin degeneracy present in the ground- 
state phases of localized electrons. In this paper we address an analogous 
problem. In our case however, the direct quantum fluctuations introduced 
by the hopping of the ions are replaced by the indirect ones, introduced by 
the hopping of an extra quantum system of free electrons coupled to the 
classical system of ions in a Falicov-Kimball way. In this new set up, the 
effect of quantum fluctuations due to the hopping could have been studied 
rigorously. In Section 3, we have considered only the hole-particle symmetry 
point. For any strength of the coupling between electrons and ions and any 
hopping, we obtained that the ground-state configurations of singly occupied 
sites (when there is no hopping of electrons), become antiferromagnetically 
ordered upon switching the hopping. This result agrees with the predictions 
of the perturbation theory for the extended Hubbard model. However, far 
away of the symmetry point, inside the domains C + (C_), where without 
hopping the ground-state configurations are spin degenerate, we also find an 
antiferromagnetic order upon switching the hopping. This result is in contrast 
with the predictions of the perturbation theory for the extended Hubbard 
model. 
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Appendix 

In order to prove Theorem 1, it is enough to notice that if Q x = n l x — r/2, 
then (8) reads: 




Rewriting Hf in terms of a nearest-neighbour bond potential (/>{ x ,y) 
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and minimizing 0/ XjJ/ \ one finds that it is a m-potential, which gives us the 
periodic ground-state configurations of 77°. The set G % of configurations which 
realize the minimum of the energy E l c (W > 0), in the terms of Q x , consists 
of configurations whose restriction to the nearest- neighbour bonds (x, y) is as 
follows: 



a) r = 2k,AI>zW, (Q x ,Q y ) = (0,0), 

b) r = 2k, 47 < zW, (Q x , Q y ) = (-k, k), 

c) r = 2k + 1, 47 > zW, (Q x , Q„) = (-1/2, 1/2), 

d) r = 2k + 1, 47 < zW, (Q x , Q y ) = (-k - 1/2, k + 1/2); 

or in terms of n l x 

a) 47 > zW, r = 2k + 1, (<, n* ) = (fc, k + 1), 

b) 47 > zW, r = 2k, (<, n* ) = (fc, fc) 

c) 47 <zW, (n x ,nl) = (r,0). 

Periodizing these bond configurations one obtains G l . 
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